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Abstract,  The  vector  wave  equation  for  an  eiectromaRnetic  field 
outside  a  perfectly  conducting  sphere  situated  in  an  inhomogeneous 
anisotropic  medium  is  reduced,  by  means  of  a  dyadic  Stratton-Chu 
formula,  to  an  equivalent  vector  integral  equation.  The  anisotropy  is 
presumed  to  arise  from  a  magnetic  dipole  at  the  center  of  the  sphere, 
so  I  hat  the  kernel  of  the  integral  equation  consists  of  the  inner 
product  of  the  appropriate  conductivity  tensor  and  the  Green's  dyadic 
(in  a  form  due  to  C.  T.  Tai)  for  a  sphere  in  free  space.  A  discussion 
is  given  of  the  spherical  vector  wave  functions  involved,  and  various 
transformations  are  applied  to  render  the  vector  integral  equation  more 
tractable.  Finally  the  vector  system  is  reduced  to  a  single  scalar 
integral  equation  apparently  more  suited  to  numerical  solution  through 
proper  redefinition  of  the  domain  of  integration. 


1.  Introduction.  The  vector  wave  equation  of  interest  has  the  form 

curl  curl  E  =  k  iK> E  (l) 

~  o  .. 

where  E  is  the  unknown  field,  k^  the  free  space  wave  number,  and 
a  dyadic,  or  second  order  tensor,  that  expresses  the  anisotropy  of 
the  medium.  Such  a  differentisd  equation  arises  naturally  in  several 
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contexts  [13,  but  in  the  present  instance  interest  is  centered  on  the 
electromagnetic  field.  Eq.  (1)  comes  about  from  Maxwell's  equations 
in  the  following  way.  Maxwell's  equations,  in  a  form  appropriate  to 
an  Inhomogeneous  anisotropic  medium,  read 


curl  E  =  p 


dH 

b 


\{f 

"  ‘o  5t  =  5? 


(2) 


£eff  =  ° 


div  B  =  0 


(5) 


where  J  is  related  to  E  by  the  conductivity  tensor 


J  =  0 


(4) 


and" the  (effective)  electric  induction  i®  given  by 


]}  ..  s  c  K  •  E 
"eff  o  ^  " 


(5) 


The  effective  dielectric  tensor  HC  is  related  to  the  conductivity 
tensor  by  the  equation 


MC  =  X  *  c  / 

o 


(6) 


Assuming  harmonic  time  dependence  substitution  of  the  second 

of  eqs.  (2)  into  the  first  then  leads  to  the  vector  wave  equation  (1) 
upon  talcing  account  of  (3). 

If  one  expands  the  left  hand  side  of  (1)  by  the  usual  vector 
identity,  a  term  in  grad  div  E  results.  In  an  inhomogeneous  anisotropic 
medium  this  term  inextricably  couples  the  field  components.  Thus  not 
only  the  vectorial  nature  of  the  boundary  conditions  but  also  the 
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Since,  for  an  arbitrary  vector  F  and  dyadic  G,  both  in 

V.  (FXVXC)=Vxr*7XG-F*\7X^xO 
and 

V  -  (<QXVXF)=CXF.7X0-7XVXF*C 

=  7XF‘7X0-0^*VX7XF, 

T 

where  G  denotes  the  transpose  of  G,  the  difference  of  the  volxine 
Integrals 

J*  V  .  (F  X  V  X  C)dV  and  J  7  *  (G  X  V  X  F)dV 
V  ^  V  ~ 

yields,  in  view  of  Gauss'  theorem,  the  following  vector  Green's  formula: 

r{F.VXVG-G’^*Vxfx  F}dV  =  fn*  (GXVXF-FX7X  G)d5:  (7) 

V  ~  8V  ~  ~ 

The  terms  in  the  surface  integral  may  be  expanded  eis  follows: 

n*GXVXFt=VXF*nXC  and  n.FxVxO  =  nxr.VXO 

and  we  finally  obtain  from  (7)  the  desired  generedlzatlon  of  the  Stratton-Chu 
formula : 

/{F.vxvxG-vx7xF*0)dV=  /{f.vxvxG-o’^'Vxvx  rJdV  = 

=^^{vxF.nXG-nXF.VX  0)dJ:  .  (8) 

The  order  of  the  inner  products  in  these  Integrands  Is  essential. 

Let  us  now  apply  formula  (8)  to  the  vector  wave  equation 
VXVxE  =  k^<IC.E  +  J(r) 

'N#  O  ^ 


f 


(9) 
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where  j(r)  represents  the  current  distribution  of  a  finite  source  that  may 
be  present, the  conduction  current  in  the  anisotropic  conducting  medium 
having  already  been  included  in  •  B  .  We  thus  take  F  =  E  and 
C  to  be  the  dyadic  Green’s  function  for  a  homogeneous  isotropic  medium; 

P  X  V<D  =  k^O  +  I  6(P  -  Q),  (10) 

but  satisfying  the  same  boundary  and  for  radiation  conditions  as  E  . 

Thus  at  the  surface  of  the  perfectly  conducting  sphere, 

n  X  E  =  0  and  n  x  G  s  0  (11) 


Substitution  of  (9)t  (10),  and  (11)  into  (8)  and  taking  account  of 
Silver-Muller  type  radiation  conditions  for  E  and  G  then  leads  to 
the  following  equivalent  integral  equations  for  the  electric  field  at 
point  P  in  an  inhomogeneous  anisotropic  medium  exterior  to  a  perfectly 
conducting  sphere; 

E  (P)  =  P  (P)  +  k^  J*  B  (Q)  •  G(P,q)dQ  (12) 

or 

E  (P)  =  F  (P)  +  k^  |c’^(P,Q)  .  aff(Q)  •  E(;i)d(4  (13) 

a 

where 


r  (P)  =  -  ;(i  <«>  • 


for  eq.  (12) 
for  eq,  (13) 


(14) 


V  denotes  that  portion  of  3~space  exterior  to  the  sphere  of  radius 

£l 


0  =  m  -  1  =o  /  (iwe^). 


a,  and 


(15) 
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It  Is  appropriate  here  to  consider  the  forms  that  (which, 

apart  from  multiplication  constants,  Is  essentially  the  conductivity 
tensor)  can  take  In  various  magnetoplasmas .  Each  form  Is  obtainable 
from  the  pondemotlve  equation  for  a  magheto-lonlc  medliim: 

3v 

®g^=q  Ce  +  vX  (B^  +  b)3  -  V  m  V  ,  (l6) 


where  q  and  m  denote  the  charge  and  mass,  respectively,  of  a  charged 
particle  moving  with  velocity  v  In  a  medium  with  ambient  magnetic 
field  and  average  collision  frequency  v.  The  form  of  conductivity 

tensor  Is  then  determined  In  the  usual  way  by  assuming  harmonic  time 
dependence,  expressing  the  vector  product  In  (16)  In  the  appropriate 
coordinate  system,  and  Introducing  the  polarization  vector  P(r)  =  (iw)  ^J(r) 
and  the  notation  of  magneto-ionic  theory  [4].  For  a  uniform  magnetic  field 
parallel  to  the  z-axls  we  have  the  familiar  result 


uniform 


X 


'  TJ  -lY  0  \ 

lY  U  0 

\o  0  (  -  Y^)/U  / 


(17) 


For  a  magnetic  vector  B^  lying  entirely  In  the  meridian  plane  In  spherical 
coordinates  the  tensor  ^  takes  the  form 


'sph*  ■  ^2  .  t2 


(18) 
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For  dipolar  coordinates  [53  P  =  cos  0/r^,  a  =  r  /  sin^©  ,  cp  =  m  , 
with  the  ambient  magnetic  field  being  directed  along  a  line, of  force 
a  =  const. (i.e.,  parallel  to  1^),  takes  essentially  the  same 
form  as  (17) : 


dipolar 


(U^  -  Y^)  /D  0 

0  D 

,  0  iY 


.(19) 


where  Y  =  (i)„  /  u. 

n  » 

We  note  that  in  all.  the  forms  (17)  -  (19)  the  conductivity  tensor 
consists  of  a  factor  X=  involving  the  electron  density, 

multiplying  a  matrix  which  expresses  only  the  effects  of  the  ambient 
magnetic  field  and  the  collision  frequency.  Thus 

X  (p)|f  ,  (ao) 

where  the  form  of  Tf  is  immediately  clear  from  the  appropriate  one  of 
expressions  (17)  -(19). 
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3.  The  Green’s  Dyadic 

A,  Tai’s  Form  of  the  Green’s  Dyadic 

The  Green's  dyadic  that  appears  in  the  integral  equations 

(12)  or  (13)  applies  to  the  exterior  problem  for  a  perfectly  conducting 

sphere  situated  in  homogeneous  isotropic  space.  The  spectrum  of  the 

linear  operator  in  (10)  is  therefore  in  general  continuous,  since  the 

volume  exterior  to  the  sphere  is  infinite.  To  determine  G  one  can 

either  employ  an  integral  representation  over  a  contour  that  guarantees 

outgoing  waves  at  infinity  or  else  adopt  a  series  expansion  in  terms  of 

vector  spherical  wave  functions,  each  term  of  which  consists  of  a 

generalized  spherical  harmonic  in  the  "angular  coordinates"  6  amd  q)  , 

multiplied  by  a  discontinuous  function  of  the  radial  coordinate.  Since 

O  in  the  latter  form  is  already  available  from  the  investigations  of 

C,  T,  Tai  [2],  it  will  be  adopted  here.  The  transpose  of  Tai’s  Green’s 

dyadic  (this  being  the  form  appropriate  to  (I3))*  properly  modified  for 
‘fiwt 

e  time  dependence,  reads  as  follows: 


where  is  the  Neumann  symbol  (1  for  m  =  0,  2  otherwise). 


c 


n 


ik 

o 


2n  1 
n(n+l) 


» 


(21a) 


(21b) 


(22) 
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and  9 

~®mn 

0 

functions 


and 
(1,4) 


mn 


He 

-fmn 

o 


and 


are  defined  in  terms  of  the  spherical  vector  wave 
as  follows: 

~  mn 
o 


0^(P) 

mn 

o 


M^<P) +  (P) 


-e 


mn 


mn 


(25) 


H  (P) 

~®mn 

o 


«J(P) 

mn 

o 


R*  (P) 
n  »e 

mn 

o 


(24) 


The  spherical  vector  wave  functions  themselves  will  be  discussed  in  the 

K  G 

next  section.  The  reflection  coefficients  R  and  R  au'e  defined 
*  n  n 

as  follows: 


a) 


S" 


In  the  sequel  we  shall  often  simply  write  r  for  the  radial  coordinate 
rp  of  P  and  r’  instead  of  r^. 

We  note  that  ^  (P,Q)  =e>«3,P),  <I%(P,Q)  =C<(Q,P),  simply 
interchanging  antecedents  and  consequents  in  the  dyadics. 


B,  Vector  Spherical  Wave  Functions. 

The  scalar  spherical  wave  functions  have  the  form  [2], [3, sec. 7, 11] 


♦^(r,n)  =  (k^r)  Pj  (cos  ©)  mcp  , 

mn 


(25) 


i 
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where  denotes  a  spherical  Bessel  function  and  0  is  the  angular 

variable  (0,cp).  We  note  that 


mn 


+ 


0 


I*  r  I* 

The  spherical  vector  wave  functions  L  ,  M  ,  N  (the  superscript 

~  mn  ~  mn  ~®mn 
o  o  o 

r  may  on  occasion  be  omitted)  are  then  defined  in  terms  of  ^  as 

®mn 

o 

follows  t 


dz  (k  r)  dP™ 

e  =  ^  ^e  =  ~dr°"  ^  ^  Ir  *  Cin^®  i©  ^ 

mn  mn 

o  o 


ffl 


r  sin  © 


z  (k  r)p”(cos  tp  1 

no  n  cos  ^  ^ 


~<P 


(26) 


=  V  X  (*  r)  =  L  X  r  = 
-e  ^e  ~  «,e  ~ 

.mn  mn  mn 

o  o 


■'  o 


=  *  ilTs  <!'  h  -  i,  <2’’) 


dP” 


~<P 


z  (k  r) 


=  k"^  X  7  X  (tg  r)  =  n(n  +  1)  °  /  pj  (cos  e){®°®}m  cp  1^  + 
~_mn  mn~  o  ~ 


[k  r  z  (k  r)]  dP® 
o  n  o  n  rcosi 

- kl -  d©"  ^sin^“  'P  i©  ^ 

o 


(28) 


[k  r  z  (k  r)]  . 

m  0  n  o  jn  /  ^xrsin-i 
+  — r— ;;  (cos  ©)l  jm  cp  1 

n  '■cos-'  ^  ^ 


sin  ©  k  r 
o 
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M  and  N  have  a  reciprocal  relation  to  each  other: 

~®mn  ~®itin 

o  o 


=  k"^  V  X  , 

-e  o  _e  ’ 

mn  mn 

o  o 


=  k"^  V  X  M*” 

«e  o  -e 

mn  mn 

o  o 


Then  we  may. write  (21)  as 


V  *  (n-m) ! 
O  (P,q)  =  2  c  2  (n+mi: 
n=l  mso 


X  (♦?*  X  (Q)  +  (Q)V']  + 


P^'^e  -Q  “^e 

mn  mn 

o  o 


+  k\x7pX  (i,f  i)«7  XV  X  Ktf  (Q)+  ^^f(Sl))r%  r'<r 
°  ^mn~  ^  ^  ^mn  ^mn 

^x  K*^(P)  +  (P))r]^X  (Q)r')  + 


♦  X  V  X  Kf  (P)  +  R®tf(t>)  )r]v  X  7  X  (tf  (q) r  0 . 

Hin  ^  ^  Mm 


r»  >  r 


This  form  of  C  (P,q)  is  helpful  in  expanding  the  dyadic  Into  components. 


The  expansion 

P|_(co6  r)  =  E  «.j°  ;  1^<C0«  9p)  V  =»  ‘S  -  ®,) 


will  play  an  important  role  in  what  follows.  It  will,  however,  be  more 
convenient  for  our  purposes  to  rewrite  (29)  in  terms  of  the  even  and  odd 
spherical  harmonics 


Y  (fl)  =  P"  (cos  e)  ^  9], 


(50) 
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Thus 


i|r  (P)  =  z  (k  Y  («), 
e  n  o  e 

mn  mn 

o  o 


(31) 


which  leads  to  an  even  and  odd  resolution  of  (29)! 


P  (cos  y)  =  P  (cos  y)  +  P  (cos  y)  = 
n  '  e,n  o,n 


(n-m)! 


m  vn+m;;  e,mn  e,mn  o,mn  o,mn 


m=o 


The  spherical  vector  wave  functions,  (30)  and  (51),  become 

bY 


bz  (k  r) 

h  «■>  =  -ir^ 

mn 

o 


-1  /,  ,  o™  ,  Oi™  , 

J  <nn^  »  r  z_^(k^r)  -55-  ig  »  -jr-  <35) 

mn  Y 


3Y 


bl 


1  e«i»,  mn 

1  omn  o  1 


!l?e  "  *n^‘'o*'^  ^sin  0  Sep  i©  3©  i(p 

mn  ^ 


(54) 


N  (P)  s  k  r  *  1  +  2k"^  V*  +  k"^(r  .  V)  V  ♦ 

~e  o  e^r  o  e  o«  e 

mn  mn  mn  mn 

o  o  o  o 


3Y 


(35) 
3Y 


(k  r)  ^{n(n+l)z  (k  r)Y  (0)1  +  [k  rz(k  r)]  [-5^ — t — S  “3 - ^3) 

o  n  o  e  -r  o  n  o  o0  _w  sin  ©  ocp  «m 

mn  Y  -T 


We  are  now  in  a  position  to  expand  the  dyadic  products  in  tl^(P,Q), 


Thus 


dY  3Y 
Smn  o""'' 


!e  <WJe  <«>  =  - 

""  Sin  ©  sin  © 


mn  mn 
o  o 


3y  By 


-1^-7 


^mn 


Smn 


9Y  3Y 

1  o™n  3mn 


^  T  T  A  w -  - 

©  Sin  ©  efp  3©’  Zep'  ”  Z(p  sin  ©'  “3©  5^  ;;©•  * 


31  3Y 

.  Qitn  Smn  , 


-cp 


~cp’ 


(36) 
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and,  defining 

y_(k  r)  =  [k  z  (k  r)] 
"^n  o  0  n  o 


,  4(v'>  =  cv W'>’' 


dY 

e 

nn 

+  n(n  ♦  i®'  ^ 

nn 

o 


ai 


ran 

‘  Ir  STS' V">  V-  * 


ai. 


nn 


+  n(n  +  l)4^V^^n^V’^^“S5^  i©  ^ 

ai 

e 

mn 

^  ^  _nn 


ai  ar 

®nn  ®mn 
o  o 


♦>h<V^^n^V'^*^io  i©’  ^ 


+  1 


ai  ai 

e  e 

nn  nn 
1  0 


=©  Bln  ©'  a© 


“5V^  V  ■" 


ai  ai 

e  e_ 


ai  aY 

e  e 

nn  nn 
o  o 


^  ■  -^-5^  .  o 

^  sin  6  Stp  a©*  ~0’  ~(p  sin  ©  sin  ©'  {p' 


(38) 
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The  summations  on  m  in  (21)  can  then  be  expressed  as 


9  P(co6  y) 
e 

mn 
o 


l!eW  !!.<«>  •  ri„  fX  is.  - 

o  0 


(39) 


a^p 


-  1  ^2^  1  .1  a®sg^  1+1  _ ^  1  ] 

>0  sin  Q  ttf*  tip  sin  ©'  i©  *  i®  aea©'  iq,"* 


and 


^  £,«>  •  (k  r)^(k  r.)P,(oo.  y)1  t 

m=o  mn  mn  n  ” 


ap 


1  •  0* 


*  “<»*»»„«‘„^o'"»r-d^is.  *  ir  nrsT  SF  ♦ 


ap 


dP 


+  n(n+l)>Jk^r)Z^(k^r')llg  "STih'*  ip.  SiT©  "^3“^  ^ 


a^P 


a^p 


^>^n^V^^n^V’^©  a©a§^  i©'  +  ^©  s©5^  ip- 


a^P 


a^p 


+  1 


1  1  .  1  J- _ O  ^  ,1 

-tp  sin  ©  atpd©*  ^©»  icp  sin  ©  sin  ©•  5^5^  ip'-*^ 


(40) 
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We  can  thus  write  (21)  as  a  dyadic  wherein  the  summation  on  m  no  longer 
appears  explicitly: 


C^(P,Q)  =  S  ^ 

<  n-=l  n 


a^p 


*  io^sin  e  sin  ©'  5^5^  ■*■  ^0^0 •  ^io»  “ 


d^P 


d^P 


1  o”  1  o” 

“  i«^sin  e  5^5©^  sin  ©'  5©5y 


d^p 


d^P 


d^P 


1  o”  1  “ 

■  iqj^sin  ©•  5©5^  "  sin  ©  5®deT  ^i©'  ■*■ 

d^P^ 
o" 

*  icp*-'^l^’^<’*‘>^d©d©»  sin  ©  sin  ©• 

•/  ®  O" 

I  ■*■  ®n^n^o**<^  H>**>^^rvr’  i©*^i-^  sln^  * 


5^ 

r,  dP^ 


+  n(n+l)(k^r^^)“^' 


dP 


dp 

8n 


^n*^o^<^  *  n^o*’>^©  dO  ip  sin^  ir^^’ 


r  ^  ^  V 


dp 


dP 

e 

n 

0 


nV>5^^o*’<^  *  5^i©'  *  Ir^  si^©'  icp'^ 


dP 


dp 


+  /V  r  Xij  <k  r  )  +  R®h  (k  r^BO-o  "v-^  1  ,  +  1  ---  ■  ^---l  , 
n  o  >^n  o  <  n  n  0  <  ~©  dO  _r^  ^  sin  ©  d^  ir^ 


\ 


r  L  <  r. 


(4lb) 
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where 


(42) 


4,<r<,r>)  .  (lcfr^;,)-Vi“(V<>  *  "X"''' V<>C'< V>>-  <'*5) 


(4)/,.  _  xTt(4), 


4,  Transformations  of  the  Vector  Integral  Equation 
A.  Reduction  to  Symmetric  Form 

It  will  now  be  assumed  that  the  function  X(P)  (i.e.,  the 
plasma  frequency)  exhibits  radial  symmetry,  and  further  that  X(r)  vanishes 
for  r  >  r^.  Then  the  Integral  equation  (13)  takes  the  form 


* 

E(P)  *  F(P)  +  f  J*  dr*  r'^  X  (r»)  J(cJ(P,<i)  •  l?(Q)  *  B((l)dfi*  ♦ 

“  ~  »  n* 

r^ 

+  /  dr'r'^  X  (r*)  JcJ(P,Q)  •  •  E(Q)dn*  ) 

**  n* 


(44) 


Evaluation  of  F(P)  according  to  (14)  for  a  horizontal  (i.e,, 
parallel  to  1^)  current  element  of  moment  p  at  the  point  (b,0^,^) 
above  the  sphere: 


J„(r,n)  =  i«p  -  b)6(e  >  <>i)6(y  - 

2  .  «  io 


r^sin  0 
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yields,  for  horizontal  polarization, 


For  vertical  polarization,  a  vertical  current  element  at  (b,Oj^,«Pj^)  , 


J„(r,n)  =  iup  -  b)6(6  -  9i)i,(T  -  Ti)  j  , 
<»v  p  ~r 

r^sin  0 


leads  in  the  same  way  to 


,  2 

K  P  CD  n  /  V  , 

=  -T 

0  n=l  m=o 


(1)  .  ,  M  (4) 

O  0 


.mn 


(46) 


(4)  Jl)  ^  J4) 

(P)  *  rK  (P)  3  ,  r  <  b 


mn 


.mn 


An  attempt  to  expand  E(P)  in  the  usual  way  in  a  series  of 


spherical  vector  wave  functions 


is  foredoomed  to  failure,  since  the  iterated  summation  causes  the  expansion, 
coefficients  in  each  row  of  the  infinite  system  of  algebraic  equations  re¬ 
sulting  from  the  substitution  of  (4?)  into  (44)  to  be  different.  Thus  unless 
some  device  (like  a  generalization  of  the  Watson  transformation  )  is  available 
that  permits  one  to  discard, for  each  q  ,  all  but  the  first  two  coefficients, 
it  will  not  be  possible  to  solve  for  the  expansion  coefficients  in  (47)  in 
the  usual  way. 


‘f 


It  may  be  possible  to  carry  out  such  an  expansion  in  terms  of  the 
spherical  Bessel  functions  and  the  Legendre  polynomials  (52)  in  y  ,  However, 
even  then  it  would  be  more  convenient  to  have  (44)  in  symmetric  form  to  assure 
the  applicability  of  the  usual  theorems  on  eigenvalues,  eigenfunctions  and 
developability  for  symmetric  kernels  [6],[7].  Thus  we  define  the  new  (definite) 
unknown  function  auid  source  function 


BiP)  =  .  E(P)  and  ^  (P)  =  |7*(P)  •  F(P) 


(48  a,b) 


and  the  symmetric  kernel 
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1 

Then  [9|P.232]  can  be  represented  in  terms  of  matrix  Lagrange 

polynomials 


a|  ^2  4  * 


(53) 


where 


L  ■f’z  A  1  A  A  y  A  A 

^  ==  i  V 


\2  \3 


^  ^21  ^23 


^  ^31  ^32 


(54) 


with 


^11.=  *1-^.  1.2.3) 


Thus 


HS  =  -  -  -  >^21 


=  Y(1  -  Y)“^  = 


x^j  =  21(1  -  =  -  Xj, 


and 


^1  =  - 


-Y^ 


-V€ 


-Vq 

Y^ 

r 

-iY. 

lY 

Y^ 

fz-  — 


-Y^ 


A-y^ 

r  © 

'  -Y  Y 
r  © 

1 

y-y| 

-lY 

r 

lY 

r 
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1  - 


Then  the  matrix  Interpolation  polynomials  are,  according  to  (5^), 


.h,r.  =  I/I,  =  I,/I  . 


Thus,  according  to  (53), 


jrrr 


} 


(55) 


(56) 


(57) 
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or 


/■ 


n 


-  gv  V. 


(58) 


where 


■ri 


1  1 

+ 


-  I 


(59) 


For  the  sake  of  brevity  we  shall  shorten  the  notation  for  the  eleaents  of 
■01^  as  follows: 


Then,  according  to  (49),  the  elements  of  ■“  ®^>(P»Q)  are  given  by 


J^2^(P»Q)  -  ^  ^r0'*^2  ^  ^  ^O0'*^2  *  ^0:p'^3^  ^ 

*  *  °W’’3  ^ 

Kjj(P,<i)  =  -Ii(p)[a*,,,j(«)  *  0*,.,^  ♦  ♦  Ii2(p)[a*^,i,2  .  ojg,,^  t  ♦ 


*  5j(P>[a  ,i|j  *  3^,5^  ♦  <L^i55  3 
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=  ’li(P)CQrr'V^^  * 


K^^CP.Q)  =  T,2(P)CQj^,f,^(Q)  + 


K22(P,Q)  =  t,2(P)[G^^,n2(Q)  ♦ 


*23(P.Q)  »  ’l2<P)Ca^r'V*^^  ^ 


KjiCp.Q)  =  n3(P)[arr.\(<i)  + 


K32<P.«i)  =  ’l3(P)[Qrr.’>2<<i>  ^• 


K^j(P,Q)  =  t,*(P)LG^^,il3(Q)  + 


+  t,3(P)[Gj^,„3  +  aj^,i,5  +  Gj^.ng] 

®r€^2  *  %'’’?]  +  ’l4(P)CGgj.'’’l  *  °e0'’*2  "*■  ■*■ 

*  V^^‘^°Jr'’'3  ^  ^^©’”5  *  ®Jq>' V 

o^-'a*  ®lp'V  *  V'’>'°te'''i  *  “es-'a  *  °«|,-’'jj  * 

*  ''6<''>“^'’'i  *  '’^■"a  *  “w-'j  ^ 

V  °^'V  *  ’i5‘”'^°L'''a  *  '’m'Ii.  *  ®^'V  * 

*  '>6“’>t®5r'’'a  *  *  “w'"?  ^ 

°^8' V  *  istPHC'ij  *  ’’L'"?  *  ®«j'V  ♦ 


*  '6<''>f®Jr'''3  *  °^''5* 
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The  complexity  of  the  elements  (6l)  casts  doubt  upon  the 
appropriateness  of  an  attempt  to  determine  the  expansion  coefficients 
in  the  usual  way.  However,  a  numerical  approach  appears  possible,  and 
we  now  direct  our  attention  towards  a  form  of  (50)  more  suited  for 
numerical  solution. 

B.  Reduction  to  a  Single  Scalar  Integrstl  Equation 

In  the  case  of  one  independent  variable  a  well-defined  method 
exists  for  reducing  a  system  of  integral  equations  to  a  single  scalar 
integral  equation  [7,  sec. 17],  C9,P.l].  The  range  of  integration  must 
be  finite,  and  the  vector  integral  equation  is  reduced  to  scalar  form  by 
simply  redefining  the  kernel  and  other  functions  involved  over  an  extended 
range  of  integration.  It  is  possible  to  apply  a  similar  method  to  (50), 
and  we  now  proceed  to  its  derivation. 

Let  the  Independent  variables  P  and  Q  in  (50)  be  respectively 

P  =  (rp,fl)  Q=(rjj,n')  ,  (62) 

and  define  new  independent  variables 

P'  =  (p,n),  Q'  =  (p',n')  (63) 

over  the  spherical  shell-like  regions  characterized  by  the  following 
inequalities : 

r  +  (j  -  l)(r^  -  a)  <  p  <  a  +  j(r^  -  a) 

(j,k  =  1,2,3) 

r*  +  (k  -  l)(r  -  a)  <  p'  <  a  +  k(r  -  a) 
o  o 


(64) 


certain  advantages  In  that  we  have  to  do  with  a  scalar  rather  than  a  vector 
Integral  equation  and  that  the  "bookkeeping"  for  machine  computation  is 
spelled  out  explicitly  in  (66)  and  (6?). 
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